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ABSTRACT 


Independent Earth orientation measurements taken by the space-geodetic techniques of 
lunar and satellite laser ranging, very long baseline interferometry, and the Global Positioning 
System have been combined using a Kalman filter. The resulting combined Earth orientation 
series, SPACE2012, consists of values and uncertainties for Universal Time, polar motion, and 
their rates that span from September 28, 1976, to April 26, 2013, at daily intervals and is 
available in versions with epochs given at either midnight or noon. The space-geodetic 
measurements used to generate SPACE2012 have then been combined with optical astrometric 
measurements to form two additional combined Earth orientation series: (1) COMB2012, 
consisting of values and uncertainties for Universal Time, polar motion, and their rates that 
span from January 20, 1962, to April 26, 2013, at daily intervals and which are also available in 
versions with epochs given at either midnight or noon; and (2) POLE2012, consisting of values 
and uncertainties for polar motion and its rate that span from January 20, 1900, to May 22, 2013, 
at 30.4375-day intervals. 



ACKNOWLEDGMENTS 


The authors would like to thank all those involved in acquiring and reducing the Earth 
orientation measurements that have been combined here to form SPACE2012, COMB2012, and 
POLE2012. This study would not have been possible without their considerable efforts. 



TABLE OF CONTENTS 


Introduction 1 

SPACE2012 2 

Data Sets Combined to Form SPACE2012 2 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 2 

Adjustments Made to Space-Geodetic Series Prior to Combination 4 

Combined EOP Series: SPACE2012 6 

COMB2012 11 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 1 1 

Adjustments Made to BIH Series Prior to Combination 11 

Combined EOP Series: COMB2012 12 

POLE2012 15 

Discussion 17 

References 19 

Acronyms and Terms 22 


v 



vi 



INTRODUCTION 


Reference series of Earth orientation parameters (EOPs) obtained by combining 
independent measurements of the Earth’s orientation are generated annually at the Jet Propulsion 
Laboratory (JPL) in support of tracking and navigation of interplanetary spacecraft. This report 
describes the generation of the most recent such combined Earth orientation series: SPACE2012, 
COMB2012, and POLE2012. Since the procedures used to generate these most recent series are 
similar to those used to generate previous such combinations, only a brief description of their 
generation is given here. Further details regarding the approach used at JPL to annually combine 
independent measurements of the Earth’s orientation can be found in Gross (1996, 2000) and 
Gross et al. (1998). 
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SPACE2012 


Data Sets Combined to Form SPACE2012 

SPACE2012 is a combination of independent space-geodetic measurements of the 
Earth’s orientation. Table 1 lists the space-geodetic series used in generating SPACE2012, 
giving their identifiers, the number of measurements from each series that were actually 
incorporated into SPACE2012, and the time interval spanned by those measurements. Note that 
the University of Texas Center for Space Research (UTCSR) satellite laser ranging (SLR) 
Universal Time (UT) values were not used in generating SPACE2012 due to problems 
associated with separating this component of the Earth’s orientation from the effects of 
unmodeled forces acting on the satellite that cause the node of its orbit to drift (see Gross et al. 
1998, p. 217 for further discussion about this point). For similar reasons, the International Laser 
Ranging Service (ILRS) satellite laser ranging length-of-day (LOD) values have also not been 
used in generating SPACE2012. 

Since it was desirable to combine only independent measurements of the Earth’s 
orientation, only those series listed in Table 1 were used, even though other space-geodetic series 
are available from other analysis centers. When more than one series determined by the same 
measurement technique was used, care was taken to make sure that the measurements themselves 
were not included more than once. In particular, polar motion measurements from the JPL 
Global Positioning System (GPS) series were only used until the start of the International Global 
Navigation Satellite Systems (GNSS) Service (IGS) combined series EOP(IGS) 95 P 01 on 
January 1, 1995; polar motion measurements from the IGS combined series EOP(IGS) 95 P 01 
were then used until the start of the IGS combined series EOP(IGS) 95 P 02 on June 30, 1996; 
polar motion measurements from the IGS combined series EOP(IGS) 95 P 02 were then used 
until the start of the accumulated IGS Solution Independent Exchange (SINEX) combined series 
EOP(IGS) 00 P 03 on February 26, 2000; and polar motion measurements from the accumulated 
IGS SINEX combined series EOP(IGS) 00 P 03 were used thereafter. Similarly, LOD 
measurements from the IGS combined series EOP(IGS) 95 P 02 were used until it ended on 
February 25, 2000 with LOD measurements from the accumulated IGS SINEX combined series 
EOP(IGS) 00 P 03 used thereafter. In addition, UT measurements from the National Oceanic 
and Atmospheric Administration (NOAA) International Radio Interferometric Surveying (IRIS) 
Intensive UT1 series were used until it ended on December 31, 1994; measurements from the 
United States Naval Observatory (USNO) National Earth Orientation Service (NEOS) Intensive 
UT1 series were then used until it ended on December 4, 2000; and measurements from the 
Goddard Spaceflight Center (GSFC) NEOS Intensive UT1 series were used thereafter. Finally, 
polar motion measurements from the UTCSR SLR series EOP(CSR) 96 L 01 were used until it 
ended on February 3, 1996 with measurements from the ILRS combined SLR series being used 
thereafter. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The Earth orientation series listed in Table 1 were first preprocessed by removing leap 
seconds from the Universal Time (UT1) values and, when necessary, by correcting the UT1 
values to be consistent with the extended definition of Greenwich Sidereal Time (GST) as 
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adopted by the International Earth Rotation and Reference Systems Service (IERS; IERS 1997, 
p. 149). Since most of the series listed in Table 1 were already consistent with the extended 
definition of GST, this correction needed to be applied to only the NOAA IRIS Intensive UT1 
series. 

Table 1 . Data Sets Combined to Form SPACE2012* 


Data Set Name 

Data Type 

Analysis 

Center 

Reference 

Data Span 

Number 

LLR (JPL13M01; VOL, UTO) 
McDonald Cluster 

LLR 

JPL 

Williams et al. (2012) 

Oct. 5, 1976, to Jan. 5,2012 

747 

OCA 

LLR 

JPL 

Williams et al. (2012) 

April 7, 1984, to April 3, 2013 

1155 

Haleakala 

LLR 

JPL 

Williams et al. (2012) 

Feb. 10, 1985, to Aug. 11, 1990 

70 

Apache Point 

LLR 

JPL 

Williams et al. (2012) 

June 4, 2006, to Aug. 28, 2012 

413 

UTCSR (CSR96L01) 
Lageos Polar Motion 

SLR 

UTCSR 

Eanes & Watkins (1996) 

Sep. 27, 1976, to Feb. 3, 1996 

2219 

ILRSA (10MAY13; Polar motion) 
ILRS Primary Combination 

SLR 

ASI-CGS 

Sciarretta et al. (2010) 

Feb. 4, 1996, to April 25,2013 

6257 

DSN (JPL97R01; T, V) 
California- Spain Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Nov. 26, 1979, to April 15, 2013 

1028 

California- Australia Cluster 

VLBI 

JPL 

Steppe et al. (1997) 

Oct. 28, 1978, to Sep. 30, 1997 

698 

NOAA (NOAA95R02) 
IRIS Intensive UT1 

VLBI 

NOAA 

Ray et al. (1995) 

April 2, 1984, to Dec. 31, 1994 

2395 

USNO (N9903) 

NEOS Intensive UT1 

VLBI 

USNO 

Eubanks et al. (1999) 

Jan. 4, 1995, to Dec. 4, 2000 

1497 

NASA/GSFC (GSF2012a) 
NEOS Intensive UT1 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2013) 

Dec. 6, 2000, to April 25, 2013 

3830 

NASA/GSFC (GSF2012a) 
Multibaseline 

VLBI 

GSFC 

NASA/GSFC VLBI Group (2013) 

Aug. 4, 1979, to April 26, 2013 

4552 

NASA/GSFC (GSFC 1122) 
Westford-Fort Davis 

VLBI 

GSFC 

Gordon et al. (1999) 

June 25, 1981, to Jan. 1, 1984 

105 

Westford-Mojave 

VLBI 

GSFC 

Gordon et al. (1999) 

March 21, 1985, to Aug. 6, 1990 

18 

GPS (21APR04; Polar motion) 
Post-processed Flinn Analysis 

GPS 

JPL 

Fleflin et al. (2004) 

June 10, 1992, to Dec. 31, 1994 

817 

GPS (IGS95P01; Polar motion) 
IGS Final Combined (ITRF93) 

GPS 

NRCan 

Kouba & Mireault ( 1997) 

Jan. 1, 1995, to June 29, 1996 

546 

GPS (IGS95P02; Polar motion) 
IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

June 30, 1996, to Feb. 25, 2000 

1337 

GPS (IGS00P03; Polar motion) 
IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to April 25, 2013 

4795 

GPS (IGS95P02; LOD) 

IGS Final Combined (IGS08) 

GPS 

CODE 

Mireault & Kouba (2000) 

Feb. 23, 1997, to Feb. 25, 2000 

1098 

GPS (IGS00P03; LOD) 

IGS SINEX Combined (IGS08) 

GPS 

IGN 

Rebischung & Garayt (2013) 

Feb. 26, 2000, to April 25, 2013 

4806 


LLR, lunar laser ranging; JPL, Jet Propulsion Laboratory; VOL, variation of latitude; UT, Universal Time; OCA, Observatoire de la Cote d’Azur; 
UTCSR, University of Texas Center for Space Research; SLR, satellite laser ranging; ILRS, International Laser Ranging Service; ASI, Agenzia Spaziale 
Italiana; CGS, Centro di Geodesia Spaziale; DSN, Deep Space Network; T, transverse; V, vertical; VLBI, very long baseline interferometry; NOAA, 
National Oceanic and Atmospheric Administration; IRIS, International Radio Interferometric Surveying; USNO, United States Naval Observatory; NEOS, 
National Earth Orientation Service; NASA, National Aeronautics and Space Administration; GSFC, Goddard Space Flight Center; GPS, Global 
Positioning System; IGS, International Global Navigation Satellite System (GNSS) Service; ITRF, International Terrestrial Reference Frame; NRCan, 
Natural Resources Canada; CODE, Center for Orbit Determination in Europe; IGN, Institut National de l'Information Geographique et Forestiere; SINEX, 
solution independent exchange; LOD, length of day. 
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Rotational variations caused by solid Earth and ocean tides were also removed from the 
UT1 values. The effect of the solid-Earth tides was removed by using the model of Yoder et al. 
(1981); the model of Kantha et al. (1998) was used to remove the effect upon UT1 of the ocean 
tides at the Mf, Mf, and Mm tidal frequencies. Since the Yoder et al. (1981) model already 
includes a contribution from the equilibrium ocean tides, only the Kantha et al. (1998) oceanic 
corrections to the Yoder et al. (1981) model were actually removed. Also note that the Kantha 
et al. (1998) model was used to remove the effect of ocean tides on only UT1, not on polar 
motion. Ocean-tide-induced polar motion variations were not removed from any of the polar 
motion observations. Finally, the only Earth orientation series listed in Table 1 that includes the 
effects of semidiurnal and diurnal ocean tides on the Earth’s orientation is the NOAA IRIS 
Intensive UT1 series. This series included these effects by adding to the released UT1 values the 
model of Herring (1993; also see Herring and Dong 1994). Hence, the same Herring (1993) 
model was used to remove them. 

On June 30, 1996, the IGS reference frame changed from the International Terrestrial 
Reference Frame ITRF93 to ITRF94; on March 1, 1998, it changed from ITRF94 to ITRF96; on 
August 1, 1999, it changed from ITRF96 to ITRF97; on February 27, 2000, it changed from 
ITRF97 to the IGS realization of ITRF97 known as IGS97; on December 2, 2001, it changed 
from IGS97 to the IGS realization of ITRF2000 known as IGS00; on January 11, 2004, it 
changed from IGS00 to a new IGS realization of ITRF2000 known as IGbOO; on November 5, 
2006 it changed from IGbOO to the IGS realization of ITRF2005 known as IGS05; on April 17, 
2011 it changed from IGS05 to the IGS realization of ITRF2008 known as IGS08 (Rebischung et 
al., 2011); and on October 7, 2012 IGS08 was replaced by the new IGS realization of ITRF2008 
known as IGb08 (Rebischung, 2012). These changes in reference frames potentially introduce 
discontinuities into the IGS combined Earth orientation series. However, the IGS Final 
combined series EOP(IGS) 95 P 02 and the accumulated IGS SINEX combined series 
EOP(IGS) 00 P 03 used here have had each of these segments aligned to the same IGS reference 
frame. Thus, to within the uncertainty in detennining the corrections required to align each 
segment, the IGS combined series EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 used here should be 
reasonably consistent with each other. They were, therefore, concatenated with one common set 
of bias-rate corrections being determined for them, as described below. 

Adjustments Made to Space-Geodetic Series Prior to Combination 

Prior to combining the series listed in Table 1 to fonn SPACE2012, series-specific 
corrections were applied for bias and rate, the stated uncertainties of the measurements were 
adjusted by multiplying them by series-specific scale factors, and outlying data points were 
deleted. Values for the bias-rate corrections and uncertainty scale factors were determined by an 
iterative procedure wherein each series was compared, in turn, to a combination of all others. In 
order to minimize interpolation error (see Gross et al. 1998, pp. 223-225), the comparison of 
each series to its reference combination was done at the epochs of the measurements of that 
series by generating its reference combination using a Kalman filter that interpolates to (and 
prints its EOP estimates at) the exact epochs of those measurements. Also, both the bias-rate 
corrections and the uncertainty scale factors for all components of a given series were 
determined simultaneously in a multivariate approach using nonlinear weighted least squares. 
Using a multivariate approach allows the correlations between the components to be taken 
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into account when determining the bias-rate corrections and uncertainty scale factors (see Gross 
et al. 1998, pp. 225). 

All the series listed in Table 1 were included in the iterative procedure. Details of the 
iterative procedure are described in Gross (1996, 2000) and Gross et al. (1998) and will not be 
repeated here. These details include: 

(1) Use of a reference series, SPACE201 1 (Ratcliff and Gross, 2013), for initial bias-rate 
alignment; 

(2) Analysis of each data type in its natural reference frame; 

(3) Clustering of the McDonald LLR stations and, separately, the DSN very long baseline 
interferometry (VLBI) stations in California, Australia, and Spain; 

(4) Initial convergence on values for the series-specific bias-rate corrections and uncertainty 
scale factors prior to the start of outlier detection and deletion; and 

(5) Final convergence on these values after detecting and deleting all data outliers. 

At the end of the iterative procedure, relative bias-rate corrections have been determined 
that make the series agree with each other in bias and rate; uncertainty scale factors have been 
determined that make the residual of each series have a reduced chi-square near one when 
differenced with a combination of all others; and outlying data points (those with residual 
values greater than four times their adjusted uncertainties) have been deleted. A total of 74 data 
points, or about 0.19% of those combined, were thus deleted. 

For GSFC NEOS Intensive and GSFC Multibaseline VLBI, the initial alignment to 
SPACE20 1 1 was performed individually for each baseline containing at least 90 data points and 
spanning more than 1 year. Baselines containing less than one year of data or fewer than 90 data 
points were aligned as a group to SPACE201 1. Subsequent steps of the iterative procedure were 
then conducted on the total intensive or multibaseline dataset as a whole. 

Starting April 17, 201 1 the IGS changed their procedure for estimating uncertainty in the 
GPS EOP values. This introduced a change in the magnitude of the polar motion X (PMX), polar 
motion Y (PMY), and LOD uncertainties of approximately a factor of 6. To account for this 
sudden change, uncertainty scale factor values after April 17, 2011 were estimated separately 
from those of earlier epochs using a trial reference series. These trial scale factors were applied 
to the original IGS GPS data which were then used in the iterative procedure described above to 
form a final reference series. The resulting scale factors shown in Table 2 for the IGS00P03 data 
set include the contribution of both the trial scale factors and the overall iteratively determined 
scale factor. 

Finally, each series was placed within a particular IERS reference frame by applying 
to it an additional bias-rate correction that is common to all the series. This additional 
correction was detennined by first combining all the series after applying to them the relative 
bias-rate corrections and uncertainty scale factors that had been detennined for them as described 
above. This intermediate combination was then compared to the IERS combined Earth 
orientation series EOP(IERS) 08 C 04 (Bizouard and Gambis, 2011) during the interval 1984— 
2012 in order to obtain the additional bias-rate correction required to make it, and therefore each 
individual series, agree in bias and rate with the IERS series. 
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The total bias-rate corrections and uncertainty scale factors determined by the procedures 
outlined above are given in Table 2. The values of the bias-rate corrections given in Table 2 are 
the sum of: (1) all the incremental corrections applied during the iterative procedure, (2) the 
corrections applied to initially align the series with each other, and (3) the additional, common 
correction applied in order to place each series within that particular IERS reference frame 
defined by the IERS combined Earth orientation series EOP(IERS) 08 C 04. The values of the 
uncertainty scale factors given in Table 2 are the products of all the incremental scale factors 
determined during the iterative procedure. The uncertainties of the bias-rate corrections given in 
Table 2 are the la standard errors of the incremental bias-rate corrections determined during the 
last iteration. There are no bias-rate entries in Table 2 for components that were either not used 
(e.g., the UTCSR SLR UT1 component) or not available (e.g., the NOAA IRIS Intensive polar 
motion components). Note that the same IERS rate correction is applied to all the data sets, 
including those such as the Westford-Mojave single baseline VLBI series, for which no relative 
rate correction could be determined. Therefore, the rate corrections given in Table 2 for those 
series for which no relative rate corrections could be detennined are simply the IERS rate 
correction, but given in the natural reference frame of that series. In these cases, uncertainties 
for the rate corrections are not given. Also note that the entries for the bias-rate corrections in 
Table 2 for the IGS combined series EOP(IGS) 95 P 02 and EOP(IGS) 00 P 03 are the same 
prior to April 17, 2011. Since these entries were initially given within the same IGS reference 
frame, they were merged, and a common bias-rate correction was detennined for them. 

Combined EOP Series: SPACE2012 

A Kalman filter was used to combine the series listed in Table 1 after the bias-rate 
corrections and uncertainty scale factors listed in Table 2 had been applied to them. The 
resulting combined Earth orientation series, SPACE2012, consists of values (Figure 1) and la 
standard errors (Figure 2) for polar motion, Universal Time, and their rates spanning September 
28, 1976, to April 26, 2013, at daily intervals; and it is available in versions for which the epochs 
are given at either midnight or noon. Leap seconds have been restored to UT 1 , and the effects 
of the solid Earth and ocean tides have been added back to the UT1 values using the same 
models for these effects that were originally used to remove them from the raw series: namely, 
the Yoder et al. (1981) model for the solid Earth tides and the Kantha et al. (1998) model for the 
ocean tides. However, semidiurnal and diurnal ocean tidal terms have not been added to and are 
therefore not included in the SPACE2012 UT1 values. 

Figure 3 is a plot of the difference between the SPACE2012 polar motion, UT1, and 
LOD values and those of the IERS combined Earth orientation series EOP(IERS) 08 C 04. 
These two series are very consistent with each other, especially after January 1, 2000 when the 
root-mean-square (rms) of their difference is only 0.053 milliarcseconds (mas) for the 
x-component of polar motion, 0.046 mas for the v-component, 0.020 milliseconds (ms) for 
UT1, and 0.018 ms for LOD. Prior to 2000, the difference between these two series exhibits 
greater variability and even some systematic behavior, particularly in the x-component of polar 
motion. This systematic behavior is due to differences in the approaches used here and by the 
IERS to correct the bias and rate of the individual series before they are combined. 
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Table 2. Adjustments to Space-Geodetic Series 2 


Data Set Name 


Bias, 

mas 


Rate, mas/yr 

Uncertainty Scale Factor 

LLR (JPL13M01 ) 


VOL 

UT0 

VOL UT0 

VOL 


UT0 

McDonald Cluster 

-0.023 ± 0.086 

0.184 ±0.088 

0.027 ±0.020 -0.032 ±0.018 

1.485 


1.362 

OCA 


-0.036 ± 0.038 

0.002 ±0.031 

0.004 ± 0.005 0.000 ± 0.004 

2.240 


1.886 

Haleakala 


-1.906 ±1.205 

-1.459 ±0.739 

-0.353 ±0.234 -0.272 ±0.155 

2.792 


2.419 

Apache Point 


-0.612 ±0.384 

-0.300 ± 0.425 

0.031 ±0.023 0.012 ±0.025 

2.003 


2.026 

DSN (JPL97R01) 


T 

V 


T 

V 

T 


V 

Califomia-Spain Cluster 

-0.822 ± 0.024 

-0.287 ± 0.058 

0.189 ±0.002 0.172 ±0.004 

2.003 


1.444 

Califomia-Australia Cluster 

0.612 ±0.018 

0.475 ± 0.053 

0.076 ± 0.008 0.006 ± 0.024 

1.436 


1.166 

NASA GSFC (1122) 

T 

V 


T 

V 

T 


V 

Westford-Fort Davis 

8.600 ±3.801 

2.971 ±5.884 

0.841 ±0.382 0.261 ±0.589 

1.260 


0.908 

Westford-Mojave 


0.416 ±1.059 

-0.034 ± 1.651 

0.002 0.001 

2.381 


0.708 

GSFC (GSF2012a) b 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

G3GcKuRi 

0.012 ±0.102 

0.136 ±0.104 

0.202 ±0.168 

-0.033 ±0.061 

0.043 ± 0.064 

-0.004 ±0.102 

1.342 

1.357 

1.909 

HrRiWfWz 

0.637 ±0.172 

-0.353 ±0.179 

0.093 ± 0.283 

0.062 ± 0.029 

-0.1 15 ±0.030 

0.006 ±0.084 

1.178 

1.249 

1.696 

MoRiWfWz 

0.153 ± 0.1 19 

-0.331 ±0.115 

-0.182 ±0.143 

-0.060 ±0.051 

-0.203 ± 0.048 

-0.050 ± 0.065 

1.412 

1.353 

0.438 

Other networks 1 

0..025 ±0.006 

0.033 ±0.005 

0.086 ± 0.009 

0.000 ± 0.000 

-0.004 ± 0.000 

-0.006 ±0.001 

1.596 

1.591 

1.666 

NOAA (95R02) 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

IRIS Intensive 



0.194 ±0.023 



0.036 ± 0.007 



1.023 

USNO (N9903) 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

NEOS Intensive 



-0.1 15 ±0.040 



0.043 ± 0.007 



1.357 

GSFC (GSF2012a) b 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

KkWz 



0.046 ± 0.025 



-0.002 ± 0.002 



1.506 

TsWz 



0.348 ± 0.068 



-0.018 ±0.004 



1. 776 

NyTsWz 



0.343 ±0.261 



-0.020 ±0.015 



1.950 

KkSvWz 



0.291 ±0.177 



-0.018 ±0.011 



1. 777 

Other networks 



-0.310 ±0.022 



0.015 ±0.002 



1. 576 

UTCSR (96L01) 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

Lageos 

-0.395 ±0.010 

0.737 ±0.009 


0.099 ± 0.004 

0.135 ±0.003 


0.758 

0.711 


ILRSA ( 1 OMA Y 1 3 ) 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

Primary Comb. 

0.032 ± 0.006 

0.004 ± 0.005 


-0.001 ±0.001 

-0.007 ± 0.000 


3.105 

3.041 


GPS (21APR04) 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

JPL Post-Flinn - 

-0.274 ±0.016 

0.142 ±0.013 


0.019 ±0.015 

-0.030 ±0.012 


1.302 

0.948 


GPS (IGS95P01) 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

Final Combined - 

-0.151 ±0.056 

0.210 ±0.054 


0.217 ±0.019 

0.313 ±0.019 


1.305 

1.025 


GPS (IGS95P02) 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

Final Combined 

0.097 ±0.005 

0.007 ± 0.004 

-0.017 ±0.011 

-0.005 ± 0.000 

-0.001 ±0.000 

0.002 ±0.001 

2.004 

1.612 

1.814 

GPS (IGS00P03) 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

PMX 

PMY 

LOD 

before 4/17/2011 

0.097 ± 0.005 

0.007 ± 0.004 

-0.017 ±0.011 

-0.005 ± 0.000 

-0.001 ±0.000 

0.002 ±0.001 

2.004 

1.612 

1.814 

after 4/17/2011 







12.478 

9.025 

10.840 


a ) Reference date for bias-rate adjustment is 1993.0. See Table 1 footnotes, mas, milliarcseconds; PMX, polar motion X; PMY, polar motion Y. 

b) GSFC VLBI Station IDs from ftp://cddis.gsfc.nasa.gov/vlbi/ivscontrol/ns-codes.txt 

G3 = NRAO 85-3; Gc = Gilmore Creek; Hr = HRAS at Fort Davis; Kk = Kokee Park; Ku = 9m at Kokee Park; Mo = 12m at Goldstone; Ny = Ny 
Alesund, Svalbard, Norway; Ri = Richmond, FL; Sv = Svetloe, Russia; Ts = Tsukuba 32m ; Wf = Westford; Wz = Wettzell. 
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UT1-UTC (ms) Polar Motion X (mas) 


A Combined Earth Orientation Series: SPACE2012 






Figure 1 . Plots of the x-component of polar motion (la), y-component of polar motion (lb), UT1-UTC (lc), and excess length-of-day (Id) as given by the 
combined Earth orientation series SPACE2012. The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Note that the 
UT1-UTC values displayed in (lc) include the tidal variations, whereas the excess length-of-day values shown in (Id) do not. 






Polar Motion X Error (mas) 


A Combined Earth Orientation Series: SPACE2012 



Figure 2. Plots of the lo formal errors of the x-component of polar motion (2a), y-component of polar motion (2b), UT1-UTC (2c), and excess length-of-day 
(2d) as given by the combined Earth orientation series SPACE2012. The inserts within panels (2a), (2b), and (2c) show that component’s post-1984 uncertainties 
on an expanded scale with the same units: milliarcseconds (mas) for polar motion, milliseconds (ms) for UT1-UTC. 
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Figure 3. Plots of the difference between the IERS combined Earth orientation series EOP(IERS) 08 C 04 and SPACE2012 formed by subtracting the SPACE2012 
values from those of the IERS series. The difference between the x-component of polar motion is shown in (3a), the difference between the j -component is 
shown in (3b), the difference between UT1-UTC is shown in (3c), and the difference between the excess length-of-day is shown in (3d). 





COMB2012 


COMB2012 extends SPACE2012 by additionally incorporating the optical astrometric 
polar motion and UT1 series that was determined at the Bureau International de l’Heure (BIH) 
from an analysis of time and latitude observations by Li (1985; also see Li and Feissel 1986). 
This BIH optical astrometric series consists of values and uncertainties for polar motion and UT1 
that span from January 5.0, 1962 to December 31.0, 1981 at 5-day intervals. 

Data Preprocessing and Treatment of Tide-Induced Rotational Variations 

The BIH optical astrometric series was first preprocessed by removing leap seconds 
from the UT1 values and by correcting the UT1 values to be consistent with the extended 
definition of GST, as adopted by the IERS (IERS 1997, p. 149). Rotational variations caused 
by solid Earth and ocean tides were also removed from the UT1 values. The same models that 
were used to remove the tidal effects from the series combined to form SPACE2012 were also 
used to remove them from the BIH series: namely, the Yoder et al. (1981) model for the solid 
Earth tides and the Kantha et al. (1998) model for the Mf, Mf, and Mm ocean tides. However, 
since the BIH UT1 measurements represent an average value over a 5-day-long observation 
window, and since 5 days is a substantial fraction of the monthly and shorter-period tides, the 
amplitudes of these solid Earth and ocean tidal tenns were attenuated prior to their removal from 
the BIH UT1 measurements. (See Gross 1996, p. 8735 and Gross et al. 1998, pp. 226-227 for 
further discussion about this point.) 

Adjustments Made to BIH Series Prior to Combination 

The preprocessed BIH optical astrometric series was combined with the space-geodetic 
series that comprise SPACE2012 after first: (1) correcting the BIH series to have the same bias, 
rate, annual terms, and semiannual terms as SPACE2012; (2) applying a constant multiplicative 
scale factor to the measurement uncertainties of the BIH series so that its residual, when 
differenced with SPACE2012, had a reduced chi-square of one; and (3) deleting those data 
points, if any, for which residual values were greater than four times their adjusted 
uncertainties. Due to software limitations associated with the need to correct the annual and 
semiannual terms of the BIH series, the above adjustments were detennined separately for each 
component of the BIH series in a univariate approach rather than simultaneously in a 
multivariate approach as was done for the series combined to form SPACE2012. The procedure 
used to determine these bias-rate and seasonal term corrections, uncertainty scale factors, and 
outlying data points has been described before (Gross 1996, pp. 8735-8738) and will not be 
repeated here. The annual and semiannual terms of the BIH series were adjusted in order to 
correct for systematic, seasonally varying effects that are known to be present in optical 
astrometric measurements. Since the values of both the BIH series and the SPACE2012 series 
are given at midnight, interpolation error (see Gross et al. 1998, pp. 223-225) is automatically 
minimized when differencing these two series for the purpose of detennining the adjustments to 
be made to the BIH series. Tables 3 and 4 give the resulting uncertainty scale factors and values 
and la standard errors of the corrections to the bias, rate, annual, and semiannual tenns thus 
determined for the BIH series. When determining these uncertainty scale factors and the 
conections to the bias, rate, and seasonal terms, no outlying data points were detected. 
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Combined EOP Series: COMB2012 


A Kalman filter was used to combine the BIH series with the adjusted space-geodetic 
series that comprise SPACE2012 after first applying to the BIH series the uncertainty scale 
factors and corrections to the bias, rate, annual, and semiannual terms that are given in Tables 3 
and 4. The resulting combined Earth orientation series, COMB2012, consists of values (Figure 4) 
and la standard errors (Figure 5) for polar motion, Universal Time, and their rates that span from 
January 20, 1962, to April 26, 2013, at daily intervals and is available in versions with epochs 
given at either midnight or noon. Leap seconds have been restored to UT1, and the effects of the 
solid Earth and ocean tides have been added back to the UT 1 values using the same models for 
these effects that were originally used to remove them: namely, the Yoder et al. (1981) model for 
the solid Earth tides and the Kantha et al. (1998) model for the ocean tides. The full amplitude 
(i.e., no tidal terms attenuated) of the effects of the solid Earth and ocean tides at the epoch of the 
time tag were added back to the UT1 values. Semidiurnal and diurnal ocean tidal terms have not 
been added to and are, therefore, not included in the COMB2012 UT1 values. 


Table 3. Adjustments to Bias, Rate, and Stated Uncertainty of Optical Astrometric Series* 


Data Set 

Bias, mas 



Rate, mas/yr 


Uncertainty Scale Factor 


PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

PMX 

PMY 

UT1 

BIH 

-2.127 ±0.860 

-0.844 ± 0.707 

-7.829 ± 1.081 

1.061 ±0.531 

1.142 ±0.433 

5.565 ±0.676 

1.828 

1.696 

1.949 

ILS 

-50.330 ±2.383 

-2.834 ± 1.691 


0.429 ± 0.487 

-0.075 ± 0.345 


2.184 

1.525 



*Reference date for bias-rate adjustment of BIH series is 1980.0. Reference date for bias-rate adjustment of ILS series is 1970.0. 


Table 4. Adjustment to Annual and Semiannual Tenns of Optical Astrometric Series* 


Data Set 


Coefficient of Sine Term, mas 


Coefficient of Cosine Term, mas 


PMX PMY UT1 


PMX PMY UT1 


BIH 

annual -5.527 ± 1.060 -7.393 ± 0.892 5.293 ± 1.394 -2.717 ± 1.132 9.403 ± 0.933 -0.908 ±1.449 

semiannual 2. 194 ± 1.083 -0.452 ± 0.907 -0.423 ± 1.421 0.981 ± 1.106 0.898 ±0.914 1.712 ±1.421 


ILS 

8.183 ± 3.361 —10.545 ± 2.381 

2.265 ±3.350 1.640 ±2.371 


annual -2.244 ±3.325 8.198 ±2.351 

semiannual 0.184 ±3.334 8.838 ±2.359 


*Reference date for adjustment of annual and semiannual terms of BIH series is 1980.0. 
Reference date for adjustment of annual and semiannual terms of ILS series is 1970.0. 
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A Combined Earth Orientation Series: COMB2012 






Figure 4. Plots of the x-component of polar motion (4a), y-component of polar motion (4b), UT1-UTC (4c), and excess length-of-day (4d) as given by the 
combined Earth orientation series, COMB2012. The discontinuous changes in the plot of UT1-UTC are caused by the presence of leap seconds. Prior to the 
introduction of leap seconds in 1972, Coordinated Universal Time (UTC) was adjusted by introducing step and rate changes designed to keep it close to UT1 
(e.g., IERS 1997, Table II-3), the effect of which is also readily apparent in (4c). Note that the UT1-UTC values displayed in (4c) include the tidal variations, 
whereas the excess length-of-day values shown in (4d) do not. 
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POLE2012 


No optical astrometric observations taken at the stations of the International Latitude 
Service (ILS) were used when creating the BIH optical astrometric series used in COMB2012 
(Li 1985; Li and Feissel 1986). The ILS polar motion measurements (Yumi and Yokoyama 
1980), which are based solely upon latitude observations made at the ILS stations are, therefore, 
independent of those comprising COMB2012 and have therefore been combined with them to 
form POLE2012. Being based solely upon latitude observations, the ILS series contains no UT1 
measurements, but consists solely of polar motion measurements that span 1899.8-1979.0 at 
monthly intervals. Although no uncertainties are given with the individual polar motion values, 
the precision with which they have been detennined is estimated to be 10-20 mas (Yumi and 
Yokoyama 1980, p. 27). An initial uncertainty of 15 mas was, therefore, assigned to each of the 
ILS polar motion values. Since this assigned measurement uncertainty will be adjusted later, its 
initial value is arbitrary as long as it is not zero and serves merely as an a priori estimate to be 
used in the series adjustment procedure described below. 

The ILS series was combined with COMB2012 to form POLE2012 after: 

(1) correcting the ILS series to have the same bias, rate, annual, and semiannual 
terms as COMB2012; 

(2) applying a constant multiplicative scale factor to the measurement 
uncertainties of the ILS series so that its residual, when differenced with 
COMB2012, had a reduced chi-square of one; and 

(3) deleting those data points, if any, whose residual values were greater than four 
times their adjusted uncertainties. 

Again, due to software limitations associated with the need to correct the annual and semiannual 
terms, these adjustments were determined separately for the x- and v-components of the ILS 
polar motion series by fitting a bias, rate, and these seasonal terms to the difference of the ILS 
series with COMB2012 during 1962.0 to 1979.0. The measurement uncertainties of the ILS 
polar motion values were adjusted by determining and applying a scale factor that made the 
residual of this fit have a reduced chi-square of one. During this procedure to determine 
uncertainty scale factors and bias, rate, and seasonal tenn corrections, no outlying ILS data 
points were deleted since no data points had residual values greater than four times their adjusted 
uncertainties. Tables 3 and 4 (in the COMB2012 section) also give the resulting uncertainty 
scale factors and values and la standard errors of the corrections to the bias, rate, annual, and 
semiannual terms thus determined for the ILS series. 

A Kalman filter was then used to combine the ILS series with the adjusted BIH and 
space-geodetic series that comprise COMB2012, after applying to the ILS series the uncertainty 
scale factors and corrections to the bias, rate, annual, and semiannual terms that are given in 
Tables 3 and 4. The resulting combined Earth orientation series, POLE2012, consists of values 
(Figure 6a and 6b) and la standard errors (Figure 6c and 6d) for polar motion and its rate that 
span from January 20, 1900, to May 22, 2013, at 30.4375-day intervals. 
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A Combined Earth Orientation Series: POLE2012 
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Figure 6. Plots of the x-component of polar motion (6a), the y-component of polar motion (6b), the Ig formal errors of the x-component of polar motion (6c), 
and the la formal errors of the y-component of polar motion (6d) as given by the combined polar motion series POLE2012. 







DISCUSSION 


The Kalman filter that was used here to combine Earth-orientation measurements 
contains a stochastic model of the process that is used to propagate the state vector and its 
associated state covariance matrix forward in time to the epoch of the next measurement. For 
polar motion excitation, the stochastic model includes a random-walk term with equal noise 
forcing both the x- and v-components of excitation and originally having a white-noise power 
spectral density of 246.4 mas 2 /day (Morabito et al. 1988). This level of polar motion excitation 
process noise in the Kalman filter was increased to 739.2 mas 2 /day in order to better match the 
observed spectrum of polar motion excitation. Increasing the excitation process noise reduces 
the level of smoothing applied to the polar motion components of the propagated state vector and 
increases the covariance associated with those components. Thus, the SPACE2012, 
COMB2012, and POLE2012 polar motion and polar-motion excitation values are not as heavily 
smoothed as were those of SPACE2003, COMB2003, POLE2003, and earlier combinations 
produced at JPL; and the uncertainties assigned to the SPACE2012, COMB2012, and POLE2012 
polar motion and polar-motion excitation values are somewhat larger. 

The LOD components of the IGS Final combined EOP(IGS) 95 P 02 and EOP (IGS) 00 P 
03 series were incorporated into JPL’s combinations starting with SPACE2007 and COMB2007. 
In general, GPS LOD estimates are contaminated by orbital artifacts that must be removed prior 
to combination with other LOD or UT1 measurements (Chin et al., 2009). When determining 
the IGS Final combined LOD series, however, the IGS applies constraints that have the effect of 
removing these orbital artifacts. Therefore the IGS Final combined LOD series was incorporated 
into the Kalman filter as a true LOD measurement whose uncertainties are assumed to be 
“white” (i.e., random, uncorrelated in time, and with a flat power spectral density). Including the 
daily IGS Final combined LOD series in JPL’s combinations helps to compensate for less 
frequent UT1 measurements and was found to improve the agreement of JPL’s combined LOD 
estimates with independent atmospheric and oceanic angular momentum values, especially at the 
highest frequencies. 

Since a Kalman filter has been used in generating SPACE2012, COMB2012, and 
POLE2012, the resulting polar motion and UT1 values are smoothed to a degree that depends 
upon both the spacing between the measurements being combined (which determines how far the 
state vector and state covariance matrix must be propagated forward in time) and the 
uncertainties that have been assigned to the measurements. Since improvements to the 
observing systems, both in the hardware and software and in the number of systems, have led to 
increasingly precise determinations of the Earth’s orientation, and since the time resolution of the 
measurements has generally increased in concert with the measurement precision, the degree of 
smoothing applied to the SPACE2012, COMB2012, and POLE2012 values is a function of 
time, with the earlier values being more heavily smoothed than the more recent values. 

Daily EOP values are reported in SPACE2012 since the NOAA IRIS and GSFC NEOS 
Intensive UT1 values are given at daily intervals, as are the GPS and ILRSA combined SLR 
values (although gaps exist in each of these data sets). However, prior to the start of these data 
sets, the measurements combined to form SPACE2012 are given less frequently; therefore, the 
Kalman filter used to combine these measurements also interpolates them in order to produce a 
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series of equally spaced values. In order to be consistent with SPACE2012, daily EOP values 
are also reported in COMB2012 even though the BIH optical astrometric series used in 
COMB2012 is given at 5-day intervals. Thus, SPACE2012, COMB2012, and POLE2012 are 
equally spaced series of smoothed, interpolated Earth-orientation parameters. 

The combined Earth-orientation series SPACE2012, COMB2012, and POLE2012 are 
available from JPL’s Geodynamics and Space Geodesy Group via anonymous ftp: 
<ftp://keof.jpl.nasa.gov/keof/combinations/2012> 
and upon request from the authors: Todd.Ratcliff@jpl.nasa.gov or Richard.Gross@jpl.nasa.gov. 
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Acronyms and Terms 


ASI 

BIH 

CGS 

CODE 

COMB 


DSN 

EOP 

G3 

Gc 

GNSS 

GPS 

GSFC 

GST 

Hr 

HRAS 

ID 

IERS 

IGN 

IGS 

ILRS 

ILS 

IRIS 

ITRF 

IVS 

JPL 

Kk 

Ku 

LLR 

LOD 

mas 

MO 


Agenzia Spaziale Italiana 

Bureau International de THeure 

Centro di Geodesia Spaziale 

Center for Orbit Detennination in Europe 

extends the SPACE series by additionally incorporating BIH optical astrometric 
measurements of polar motion and UT 1 

Deep Space Network 

Earth orientation parameters 

IVS station code for 85-3 antenna at NRAO Green Bank, WV 

IVS station code for Gilmore Creek, Fairbanks, AK 

Global Navigation Satellite System 

Global Positioning System 

Goddard Space Flight Center 

Greenwich Sidereal Time 

IVS station code for HRAS 

Harvard Radio Astronomy Station at Ft. Davis, TX 
Identity 

International Earth Rotation and Reference Systems Service 

Institut National de lTnformation Geographique et Forestiere 

International Global Navigation Satellite Systems (GNSS) Service 

International Laser Ranging Service 

International Latitude Service 

International Radio Interferometric Surveying 

International Terrestrial Reference Frame 

International VLBI Service for Geodesy & Astrometry 

Jet Propulsion Laboratory 

IVS station code for Kokee Park, Kauai, HI 

IVS station code for 9-m antenna at Kokee Park, Kauai, HI 

lunar laser ranging 
length-of-day 

milliarcsecond 

IVS station code for 12-m antenna at Goldstone, CA 
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NASA 

NEOS 

NO A A 

NRAO 

NRCan 

Ny 

OCA 

PMX 

PMY 

POLE 


Ri 

nns 

SINEX 

SLR 

SPACE 

Sv 

T 

Ts 

USNO 

UT 

UTO 

UT1 

UTC 

UTCSR 

V 

VLB I 
VOL 

Wf 

Wz 


National Aeronautics and Space Administration 
National Earth Orientation Service 
National Oceanic and Atmospheric Administration 
National Radio Astronomy Observatory 
Natural Resources Canada 

IVS station code for Ny Alesund, Svalbard, Norway 

Observatoire de la Cote d’Azur 

polar motion X 
polar motion Y 

extends the COMB series by additionally incorporating ILS optical astrometric 
measurements of polar motion 

IVS station code for Richmond, FL 
root mean square 

Solution Independent Exchange 
satellite laser ranging 

a combination of independent space-geodetic measurements of the Earth's 
orientation 

IVS station code for Svetloe, Russia 

transverse component of Earth orientation from single-baseline VLBI 
IVS station code for 32-m antenna at Tsukuba, Japan 

United States Naval Observatory 
Universal Time 

Universal Time as determined at a particular observatory 

The principle form of Universal Time; mean solar time at 0° longitude obtained 

by correcting UTO for observatory location and polar motion. 

Coordinated Universal Time 

University of Texas Center for Space Research 

vertical component of Earth orientation from single-baseline VLBI 
very long baseline interferometry 
variation of latitude 

IVS station code for Westford, MA 
IVS station code for Wettzell, Germany 
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